Strong coupling approach in dynamical mean-field theory for strongly correlated 

electron systems 



Ihor V. Stasyuk and Andrij M. Shvaika 

Institute for Condensed Matter Physics 
of the National Academy of Sciences of Ukraine, 
1 Svientsitskii Str., UA-79011 Lviv, Ukraine 



Abstract 

We review two analytical approaches in Dynamical 
Mean-Field Theory (DMFT) based on a perturbation 
theory expansion over the electron hopping to and from 
the self consistent environment. In the first approach 
the effective single impurity Anderson model (SIAM) is 
formulated in terms of the auxiliary Fermi-fields and 
the projection (irreducible Green's function) technique 
is used for its solution. A system of the DMFT equa- 
tions is obtained that includes as simple specific cases 
a number of known approximations (Hubbard- III, A A, 
MA A, . . .). The second approach is based on the dia- 
grammatic technique (Wick's theorem) for Hubbard op- 
erators that allows to construct a thermodynamically 
consistent theory when SIAM exactly splits into four 
components (subspaces): two Fermi liquid and two non- 
Fermi liquid. The results for the density of states, con- 
centration dependences of the band energies, chemical 
potential and magnetic order parameters are presented 
for different self-consistent approximations (AA, strong 
coupling Hartree-Fock and further). 



1 Introduction 

Many unconventional properties (e.g., metal-insulator 
transition, electronic (anti)ferromagnetism) of the nar- 
row-band systems (transition metals and their com- 
pounds, some organic systems, high-T c superconductors, 
etc.) can be explained only by the proper treatment 
of the strong local electron correlations. The simplest 
models allowing for the electron correlations are a single- 
band Hubbard model with on-site repulsion U and hop- 
ping energy t and its strong-coupling limit (U ^> t): t—J 
model. Recent studies of the Hubbard-type models con- 
nected mainly with the theory of high-T c superconduc- 
tivity and performed in the weak- (U < At) and strong- 
(U ^> t) coupling limits, elucidate some important fea- 
tures of these models [Q. But still a lot of problems 
remains, especially for the U ^> t case where there are 
no rigorous approaches. 

Despite the relative simplicity of the models used for 



their description the theory of electron spectrum and 
thermodynamic properties of such systems is far from 
its final completion. The use of localized (atomic) basis 
of electron states is the general feature of the models. 
Corresponding Hamiltonians 



H 



E*< 



EE**. 

(ij) ° 



(1) 



include, on the one hand the electron transfer (hopping) 
Uj between neighbouring sites (atoms) in the crystal lat- 
tice and on the other hand the short-range single-site 
electron correlations. It is primarily the on-site energy of 
Coulomb repulsion U in the case of the Hubbard model 
and the models based on that one: 



Hi = Urnimi - j»E j 



(2) 



Models like (QJ, @ can be solved exactly in two lim- 
iting cases: atomic limit (t = 0) and band electrons 
(U = 0). Near these extreme cases the expansions in 
terms of t or U are used, but the consistent formulation 
of the perturbation theory especially in the case of strong 
coupling is not a simple task. The case of an interme- 
diate coupling t ~ U is more complicated for consider- 
ation. In this region of parameter values, the splitting 
in the band electron spectrum and the metal-insulator 
transition takes place. 

Due to the presence of strong electron correlations, the 
state of the electron system and its properties depend 
essentially on the mean electron concentration. At a 
different filling of electron states and depending on the 
relation between t and U parameters, the system can 
be paramagnetic or the transition into ferro- (antiferro- 
) magnetic phase can take place. In the case of the 
more complicated structure of the Hamiltonian Hi (due 
to the allowance for the other, besides electron, degrees 
of freedom) or when the interaction is extended to the 
nearest neighbours in a lattice, the charge ordering can 
appear; the effects of phase separation become possible 
as well. The listed phenomena are the subject of study 
in the framework of various approaches and methods. 

A new impulse in the investigations in this field is con- 
nected with the development of a new approach having 
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its origin in works OEIQ] where the study of the 
©-type model in the limit of infinite dimensionality of 
space (d = oo) has been proposed. Due to the princi- 
pal simplifications in the perturbation series taking place 
in this case the possibility exists to obtain exact results 
using the scheme that corresponds to the well known 
coherent potential approximation (CPA) in the theory 
of disordered crystalline alloys. The rapidly developing 
corresponding method became known as the dynamical 
mean-field theory (DMFT). 

The central point in this method is formulation and 
solution of the auxiliary single-site problem. An initial 
model is mapped on that one while considering one site 
characteristics of the electron spectrum, such as single- 
site electron Green's function (see [UEKI, as well as the 
reviews [HUHI)- In this case the separated lattice site 
is considered as placed in some effective environment. 
Since the processes of electron hopping from the atom 
and returning into the atom are taken into account, the 
mean field acting on the electron states of the atom pos- 
sesses a dynamical nature. This field is described by the 
coherent potential J a {uj) that should be determined in 
a self-consistent way. The analytical properties of the 
solutions in the DMFT are considered in 

Only in some simple cases the single-site problem can 
be solved analytically |HU12j (e.g., the Falicov-Kimball 
(FK) model QB]). In general, including the Hubbard 
model, the application of numerical or seminumerical 
(such as quantum Monte Carlo [fi HTHITfi] or exact diag- 
onalization 1 1 fill 1 Tj as well as numerical renormalization 
group [lcf], see P]) methods turns out to be necessary. 

At present, the Dynamical Mean-Field Theory is ap- 
plied to investigate different effects in the various sys- 
tems described by the simplified or realistic models. 
First applications were devoted to the investigation of 
the changes of the density of states at the metal-insulator 
transition and appearance of the antiferromagnetic and 
ferromagnetic states in the Hubbard model [6l ll5lll9ll2f)l 
[2T1f22lf23l[2lf25ll2fill27lf28l[29] and other strongly cor- 
related electron models: Hubbard model with orbiral 
degeneracy (2DJ an d disorder |31| . boson-fermion model 
|32j . extended Hubbard model p^l, double-exchange 
model pi], Falicov-Kimball model with correlated hop- 
ping jSSj, two-band Hubbard model [36137]. periodic An- 
derson model [38| . 

Besides, different types of the response functions and 
transport coefficients are also calculated. The gen- 
eral basics how to derive response functions in DMFT 
are given in |1 111391112] . Investigations of the charge 
and magnetic susceptibilities revealed also chess-board 
charge-density-wave phase at half-filling pi] as well as 
incommensurate order and phase separation at other fill- 
ings |4()U41j . DMFT is used also to investigate optical 
conductivity |2QII21ll4^1l43j . electronic Raman scatter- 
ing pHlllfi] . thermoelectric response |4*fi] . 



In last years, DMFT is used as an approximation 
scheme to consider the electron-electron interaction to- 
gether with band degeneracy and lattice structure of the 
actual materials within the so-called LDA+DMFT ap- 
proach, that allows to describe correctly the insulating 
state of the transition- metal oxides, band structure and 
phase diagrams of the different compounds @ZEE|- 

At the same time it is of interest to develop approx- 
imate analytic approaches to the solution of the single- 
site problem. Their application at that stage is more 
effective than at considering the full model (a short re- 
view of such attempts was given recently in |51U52j l. 
The availability of the analytical (even of approximate) 
method is useful especially for new models as well as at 
the transition to the finite dimensionality of the system. 
The accuracy of approximation can be estimated relat- 
ing to the results of numerical calculations. 

The first analytical approximation proposed for the 
Hubbard model was a simple Hubbard-I approximation 
|53] (see Ref. p2] for its possible improvement) which 
is correct in the atomic (t — 0) and band (U = 0) 
limits but is inconsistent in the intermediate cases and 
cannot describe the metal-insulator transition. Hub- 
bard's alloy-analogy solution [&£] (so-called Hubbard-III 
approximation) incorporates into the theory an electron 
scattering on the charge and spin fluctuations that al- 
lows us to give qualitative description of the changes 
of the density-of-state at the metal-insulator transition 
point. Hubbard-I and Hubbard-III approximations in- 
troduce two types of particles (electrons moving between 
empty sites and electrons moving between sites occu- 
pied by electrons of opposite spin) with the different en- 
ergies that differ by U and form two Hubbard bands. 
Related schemes of the so-called two-pole approxima- 
tions jUllHZj, which are justified by the t/U -C 1 per- 
turbation theory expansions [SHI ^ are a ^ so considered. 
However, in the recent QMC studies (SUED] there are 
clearly distinguished four bands in the spectral func- 
tions rather than the two bands predicted by the two- 
pole approximations. Such four-band structure is re- 
produced by the strong-coupling expansion for the Hub- 
bard model jnOl m the one-dimensional case. There are 
also analytical approximations developed specially for 
the effective environment in the DMFT |4"5ICT| . Within 
other approaches let us mention non-crossing approxi- 
mation j2l]|S], Edwards-Hertz approach OEj], iter- 
ative perturbation theory jEUESj, alloy-analogy based 
approaches ;(iO 07 , and linked cluster expansions [4llfl8| , 
which are reliable in certain limits and the construction 
of the thermodynamically consistent theory still remains 
open |52] , 

The aim of this paper is to review two recently pro- 
posed approaches 69, 70^ based on the rigorous pertur- 
bation theory scheme in terms of electron hopping for 
the Hubbard-type models. 
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The first approach is based on the technique of 
the irreducible Green's functions. The procedure of pro- 
jecting onto the basic set of operators is used (the set 
consists of the single-site electron Hubbard operators of 
the Fermi-type). The recipe is given for the construc- 
tion of the system of equations for the coherent potential 
and self-consistency parameter (having the meaning of a 
static part of the effective internal field) in the approach 
that is a generalization of the Hubbard-III approxima- 
tion. Specific cases are considered corresponding to the 
more simple approximations of the alloy-analogy (AA) 
or modified alloy- analogy (MA A) type OEE] in the 
DMFT method as well as to the certain decoupling pro- 
cedure in the two-time Green's function method when 
applied to the initial electron problem. 

Another possibility is to build analytical approaches 
by the systematic perturbation expansion in terms of the 
electron hopping using diagrammatic tech- 

nique for Hubbard operators [731175] . One of them was 
proposed for the Hubbard (U = oo limit) and t— J mod- 
els |76U77j . The lack of such approach is connected with 
the concept of a "hierarchy" system for Hubbard opera- 
tors when the form of the diagrammatic series and final 
results strongly depend on the system of the pairing pri- 
ority for Hubbard operators. On the other hand it is 
difficult to generalize it on the case of the arbitrary U . 

In the second part of this paper we show how a rig- 
orous perturbation theory scheme in terms of electron 
hopping that is based on the Wick's theorem for Hub- 
bard operators |74U75j and is valid for arbitrary value 
of U (U < oo) and does not depend on the "hierar- 
chy" system for X operators can be developed for the 
Hubbard-type models (Tj]] . In the limit of infinite spatial 
dimensions, these analytical schemes allow us to build 
a self-consistent Baym-Kadanoff type theory |78lf7 9j for 
the Hubbard model and some analytical results are given 
for simple approximations. The Falicov-Kimball model 
is also considered as an exactly solvable limit of Hubbard 
model. 

2 Hubbard model and similar models in 
a limit of infinite dimension of space 

(d = oo) 

The transition to the d = oo limit in the DMFT ap- 
proach is accompanied by the scaling of the electron 
transfer parameter 

t* 

t= V3' (3) 

In the case of d-dimensional hypercubic lattice with an 
electron spectrum 



Sk 



d 



2f 



V Oi = l 



cos k a a , 



(4) 



this procedure leads to the Gaussian density of electron 
states |2] 



i 



■ exp 



to 
"4^2 



(5) 



The average kinetic energy remains constant in this case 
in the limit d = oo. 

The scaling © has a significant effect on the struc- 
ture of diagrammatic series for single-electron Green's 
functions of the model of the and J2J type. In partic- 
ular, the irreducible self-energy part of such a function 
becomes a purely local (a single-site) quantity 01 



Stj,tr(w) = £ CT (w)<^ 



(6) 



The Fourier-transform of £ij )CT (w) is hence momentum- 
independent 

S <r (k,w)=S <7 (w). (7) 

This leads to tremendous simplifications in all many- 
body calculations for the Hubbard model and related 
models and enables us to obtain the exact numerical 
results for the main parameters of the electron spectrum, 
to describe magnetic phase transitions and the metal- 
insulator transformation etc. (see, for example, [9ll52]l. 

The possibility of obtaining exact solutions in the d = 
oo limit opens the way to the development of a theory 
based on the expansion in powers of 1/d (the results for 
d = oo can be considered as zero approximation in this 
case). Such approaches have been elaborated for the 
last few years [80.81J. On the other hand, consideration 
in the framework in the d = oo limit is not only of an 
academic interest. It turns out that a set of the known 
approximating schemes or methods is correct in the d = 
oo limit. Besides, the obtained physical conclusions can 
be transferred in many cases to the system with finite 
dimensions keeping their suitability even at d = 3. 

The formal scheme of calculating the electron Green's 
functions and the main thermodynamical quantities can 
be developed basing on the diagrammatic expansions in 
powers of interaction parameters (such as energy U in 
the case of the Hubbard model) or matrix elements of 
the electron transfer tij. The electron Green's function 
in (k, uS) representation 



(8) 



can be expressed in the first or in the second of these 
cases as 

G kH = , , „ r ^ 7, * (9) 



G£M = — 



S CT (uj) - t k 



(10) 
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where S CT (w) or E a (uj) are the irreducible paxts (in 
the diagrammatic representation) according to Dyson or 
Larkin, respectively, 



S G . 1 (w) = 10 + (i - Y, a {u) . 



(11) 



To calculate the Y, a {ui) [or B, a (u>)} function, the effec- 
tive single-site problem is used. As was shown in 
the transition to this problem corresponds to the replace- 
ment 



(12) 

xTexp |- jf dT jf dr' ^ J ct (t-t')4«MO j , 
where 

Ho = Hi (13) 

and J a {r — r') is an effective auxiliary field which is 
determined self-consistently from the condition that the 
same irreducible part E a (uj) determines the lattice func- 
tion tff]~0fi as well as the Green's function (w) of the 
effective single-site problem. The last one is connected 
with 2o-(w) and J a (u>) by the relation 

GV>{w)= 1 (14) 

On the other hand, 

GW( W ) = G iijCT M = i^G£(u;). 



(15) 



Dynamical field J a {r — t') describes electron hopping 
from the given site into the environment and vice versa; 
the electron propagates in the environment without go- 
ing through this site between moments r and r' . The 
expression 

J.M = 5> fe %Gg CT (16) 
kj 

corresponds to this situation (the relation 111611 is known 
from the standard CPA scheme |82I83| 1: here Gj£j g is the 
electron Green's function for a crystal with the removed 
site i. 

The set of equations (1101 . 11411 and 11511 becomes 
closed when it is supplemented by the functional de- 
pendence 

G« fl) H = /([4HD, (17) 

which is obtained as the result of solving the effec- 
tive single-site problem with the statistical operator 
exp(— /3H c ff). It is possible to do this in an analytical 
way only in some cases of simple models (a Falicov- 
Kimball model [H]; a pseudospin-electron model at 
U = [Bl|; a usual binary alloy model). In general, 
numerical methods are used. 

The scheme described lies at the basis of the above 
mentioned DMFT approach used in the last years in 
considering strongly correlated electron systems. 



3 Electron Green's functions of the ef- 
fective single-site problem 

As it was mentioned, the central point in the DMFT ap- 
proach is the solution of the effective single-site problem 
and the determination of the connection between the dy- 
namical mean field (coherent potential) J a {uj) and the 
single-site electron Green's function Ga(u>). Recently 
an approximate scheme [01 , which is based on the tech- 
nique of the irreducible two-time temperature Green's 
functions and leads to the results having an interpolating 
character, was proposed. The Hubbard model is taken 
into consideration to illustrate the method. 

Let us reformulate a single-site problem introducing 
explicitly an effective Hamiltonian 



H cS = H a + VJ2(aUa + &a a ) + 



(18) 



where the auxiliary Fermi-field (£<r,££) is brought in. 
It describes the environment of the selected site and 
formally is characterized by the Hamiltonian H^. The 
single-electron transitions between the site and the en- 
vironment are taken into account. 

An explicit form of the Hamiltonian H^ is unknown. 
Let us consider, however, the Green's function 



g*{u) = {{&{&)) 



(19) 



for auxiliary fermions as the given function. The Green's 
function Q a (j—T r ) = (T T £J. (t)£o-(t'))£ [where averaging 
is performed with the part iff of the Hamiltonian II18II ] 
corresponds to the function ifTfljl in the Matsubara's rep- 
resentation. It is shown in [6^ that the expansion of 
the exp(— (3H e ff) operator in powers of V and the sub- 
sequent averaging over the states of ^-subsystem using 
the Wick's theorem and functions lTl"9l) leads to the sta- 
tistical operator l(T2|: 

(expt-^fles )}<*«> = cxp(-/3ii cff ) . (20) 

The relation 

2nV 2 g a {uo) = Mlo) (21) 

takes place in this case. 

The obtained result points out to the possibility of the 
Green's function G„ (lu) calculation based on the Hamil- 
tonian H c g. The averaging over the a, a^-variables is 
performed with the use of the Gibbs distribution while 
over the £, ^-variables it is done with the help of func- 
tion dSJ. 

Let us write the Hamiltonian (TTsl) for the case of the 
Hubbard model in terms of Hubbard operators 

H cS = -J^X™ + 2X 22 \ +UX 22 + 

+VY, [(X^ + oX 2s )i a + £{X°° + aX s2 )\ + 

a 

+H 6 . (22) 
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Here the basis of single-site states \n^n{) 

|0) = |0,0), ||} = |0 ) 1) 
|2) = |1,1), |T) = |i,o) 



A 2S 



(23) 



is used (a =|, 4). In this case the Green's function 
g\j\lo) can be written in the form 

G (a) = ((X°-\X M )) u +a((X°-\X 2 -)) u 

+ a{{X s2 \X^)) u + {{X s2 \X 2s ))^ (24) 

(a representation in terms of the two-time Green's func- 
tions is used). 

We will write the equations for functions ll24^l using 
the equations of motion for A-operators: 

i±X°°(t) = [X -, H cS ] = -liX°° + V(X Q0 + x™)t; a 



dt 

d 



VX aa ^ + aVX 02 ^ 



^ t X s2 {t) = [X* 2 ,H cfi ] = (U- n)X s2 



(25) 



+ aV(X 22 + X^^-oVX^is-VX 02 ^ 

In the Green's functions of higher order we shall sep- 
arate the irreducible parts using the method developed 
in |85l8fij . Proceeding from the equations of motion ll25|i 
we express derivatives idA 0<T ( CT2 )/dt as a sum of regular 
(projected on the subspace formed by operators A n<T , 
X a2 ) and irregular parts. The latter ones describe an 
inelastic quasiparticle scattering. We obtain 



[A UCT , H cB ] = -fiX {>a + a^A u<T + a° 2 a X s2 



7 0tT 



[X s2 , iT eff ] = (U- ii)X s2 + a\ 2 X l 

'2 



_ GJ"2 v<r2 i ryrrl 



(26) 



Operators Z® a and Z° 2 are defined as orthogonal ones 
to operators from the basic subspace: 



({Z°^ s2 \X aO }}=0 , 
({Z a ^ 92 \X 2rf }) =0 . 



(27) 



0ct(ct2) 



These equations determine the coefficients a i 

Using the described procedure we come to the expres- 
sions 



Z 0a = V(X 00 + X™)£ a + VX S °& + aVX a2 tl , 
Z s2 = aV(X 22 + X ss )£ a - aVX^ts 



VX 02 £ 



(28) 



where 



(X 22 +X™)£ a = (X 22 + X 9 °)£, a 
1 



1 



X 02£T =X ^1- -j- {X° 2 d)X°° 
A()a 



1 



-^-(£x m )x* 2 

A2a 

and A pq = (Xpp + Xii)- A 0a = 1 - n s , A s2 = n 6 
Here 

Oct 5-2 V 

a l = -aa 1 = - — <p a , 

AQa 

= -GO?% = J— at P<J > 

Ala 

ip a = {&X™) + a{X" 2 £) . 



(29) 



(30) 



(We put (fia = ip*). 

The equations for the first two functions in (I24T) have 
in this case the form 



Ul — 

v 



«,-b„) {((X* 2 \X° )) 
{{Z s2 \X a0 )) 



(31) 



where the notations 



V V 
a CT = —//+- — tp<j , b a = U-fi+- — ipcr (32) 

AQa Ala 

are used (a similar set of equations can be written for 
functions ((Z 0a ^ \X 2S ))). 

An equation for the Green's function ((Z 0a ^ |A CT °)) 
(as well as for the function {{Z aa ^\X 2s ))) can be ob- 
tained by means of the differentiation with respect to 
the second time argument. Applying the similar proce- 
dure of separation of the irregular parts we obtain the 
expression 

G = Go + GoPaGo , (33) 
where the matrix Green's function 

r-^ ({{x °\x°°)){{x °\x^))\ 

U ~ Z7T ^((A^A* )) {{X S2 \X 2S )) ) [M) 



is introduced. G is a nonperturbed Green's function 

u; • ^ 



r, _ 1 ( u-ba /4 

^0 — tt" I _ V 



where 

D a = {uj — a a ){uj - b a ) 

and 

P„ = 2n 



V 2 



<£ (36) 



A^a-Aia 

\ f((Z°°\Z°°))((Z a °\Z 2 *)) 



A£) U(Z™\Z° )) ((Z™\Z 2 *)) 



'A^l 



(37) 
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has the meaning of a scattering matrix. Being expressed 
in terms of irreducible Green's functions it contains the 
scattering corrections of the second and higher order in 
powers of V. The separation in P of the irreducible, with 
respect to V, parts enables us to obtain a mass operator 
M 

P a = M a + M a G M a + M„G MaGoM* + ■■■, 
M lT = P a \ il . (38) 

In this case the relation l(33jl can be transformed into 
the Dyson equation 

G = Go + G a M<jG (39) 

with the solution 

G = (1 — G M a )^ 1 G , (40) 
which provides a final expression for the Green's function 

4 Different-time decoupling of irreduc- 
ible Green's functions 

We will restrict ourselves hereafter to the simple approx- 
imation in calculating the mass operator P, taking into 
account the scattering processes of the second order in 
V. In this case 

Ma = , (41) 

where the irreducible Green's functions are calculated 
without allowance for correlation between electron tran- 
sitions on the given site and environment. It corresponds 
to the procedure of the different-time decoupling |8?| . 
which means in our case an independent averaging of 
the products of X and £ operators. 

Let us illustrate this approximation with some exam- 
ples. 

1. The Green's function 



According to the spectral theorem we have 



+ 00 

i r du' 



1) 



At 



-iw't/tt 



2tt J lo — u>' J 2ir 

— oo 

x (x m + x aa ) t {x m + x aa )^y T 



(£(*) 



(42) 



Due to the different-time decoupling 

(£(t)(x 00 + x™) t (x 00 + x™)^y* 

« ((X m + X™) t (X m + X™))(£m«) .(43) 

We will take the first of these correlators in a zero ap- 
proximation 



((X 00 + X a<r ) t (X 00 + X™)) 



((X c 



x c 



A 



On 



(44) 



and substitution of l(i3jl into lH2)) leads in this case to 
the result 



2. The Green's function {{X Sa £, < y\^X aS )) bJ = I 2 {u). 

The representation of the I%{uS) function in the form 
analogous to l(4*2l) leads to the time correlation function 
{(,l{t)X ad \t)X 9 ° '£ 5 ) ir that can be approximated as 

(£t {t)X' 9 (t)X™&Y* a (X^ {t)X s °) (t)&) 

«<A™)<4(t)fc>. (46) 

In this case 



3. The Green's function ((X 02 ^ S \^X 20 ))^ = I 3 (w). 
The corresponding time correlation function is decou- 
pled as 

«exp[i(C/-2 M )t](X 22 )(6(t)^) . (48) 
Using this expression we obtain 



I 3 (0j) = ^(X m +X^)((4\^))^ 2! , , 



(X m -X 22 ' 



1 



duo' 



2tt J u>+2^-U-uj' 



xtanh^[-29f«ci|fo)) w /. 



(49) 



Let us mention that at the half-filling of electron states 
(when n = 1, (X 00 ) = (X 22 )) 



h{u) = -{X 22 ){{tM))u-2^ 
(X 22 ) 



(50) 



Following the described procedure and taking into ac- 
count the relation (12 H we will come to the following 
expressions for irreducible Green's functions: 

«Z 0CT K°))c= AoMlj) - R a (u) , 
{{Z s2 \Z 2s ))^=A 2s J a (u)-R a (oj) , 
((Z a °\Z 2s )) u = ((Z* 2 \Z°°)) U = R a (oj) , (51) 

where 

R a (u) = i(X 00 + X 22 )J S {U -2n-u)- (X™)Mu) 



(x oo _ X 22 } T^' SIM-"/ i0 + ) N , 



7T w + 2/i — U — UJ 



- tanh -— (52) 
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5 Basic set of equations 

Using the results obtained in the previous section we 
can write the expressions for mass operator components 
Mo-^afl. On the basis of relation 

^V) = (E G ^Kf 1+J -H' ( 53 ) 

a/3 

[which follows from QH ] and formula ijlflj) it is possible 
to determine the single-site self-energy part. We obtain 

E a (uj)= lo - e + A Qa - e-A 2s - VL a (uj) 

x (lo — e+)(bj — £_) — ujQ, a {w) 

1 ~ 1 

+ (e+A 2s + e-A 0ry )n a (uj)\ , (54) 
where e+ = U — fi, e_ = —fi and 

A 0a A 2S 

n a {u) = J a {u)-^L. (55) 



Ao a A 2s 

It should be mentioned that formula 1)5411 can be also 
represented in the form 



3"» 



A 0o 



Ao 



_L0 — £- — fl a (uj) LO — £ + — Q, cr (u;)_ 

+ & a (Lo). (56) 

The relation JHJ together with (JTHJl, H3J) and JUJ cre- 
ates a set of equations for the coherent potential J a {^), 
self-energy part S CT (a;) and Green's functions Gu itT {uj) 
and G£(w). 

It should be noted that the parameter ip a , which is 
expressed in terms of average values of the products of 
X and £ operators (formula (JHSJ), is a functional of the 
potential J a {uj). According to the spectral theorem 



V{X a0 ^^)=i 



(57) 



On the other hand, using the linearized equation of 
motion ll26)) and neglecting the irreducible parts, we can 
obtain the following set of equations 



V 



\ A 0rr 

+v 2 *f^mx 2 n) = ^j < ,, 

A 25 27T 

V*af^{UX°°)) 

A 0a 



(58) 



V 



+V(&\X*')) [lo-s + - — Va = o^J L 



A 2S 



A 2S 



It follows herefrom in the U — > oo limit 

+ oo 



V Va = -V{X^ 9 ) 



(IlO 



1 



2?r / + 1 



s A os J s (lo) 



-23 



LO — £- 



V- 



j+iO+ 



or, in the Matsubara's representation 



\LO v - £_ - V-¥? 



(59) 



(60) 



Thus, we obtain a self-consistent equation for the pa- 
rameter Lp a . 

6 Some specific cases 

Equations obtained in the previous section form an 
approximate analytical scheme of calculating both the 
single-site and the full electron Green's function in the 
framework of DMFT. Let us compare it with the stan- 
dard approximations known from literature which are 
based on the assumption of the single-site structure of 
the electron self-energy. For this purpose we will con- 
sider some specific cases. 

6.1 Hubbard-I approximation (J CT = 0, R a = 0, 
<p a = 0; Cl a = 0) 

It is the simplest approximation; renormalization of en- 
ergies of atomic electron transitions is absent and the 
scattering processes via coherent potential are not taken 
into account. The expression for the single-site self- 
energy part 



3„H 



A 0a 



Ao 



LO — S+ 



(61) 



corresponds to the Hubbard-I approximation Elec- 
tron energy spectrum described by Green's function 
G1{lo) consists in this case of two Hubbard subbands 
divided by a gap existing at any relationship between 
the values of U and t parameters. 

6.2 Static mean-field approximation (J CT = 0, 
R a = 0; Cl a = B = V<p a /A Q(T A 2S ) 

In this case only a self-consistent shift of the energy lev- 
els of the single-site atomic problem is taken into ac- 
count. The coherent potential J a (cj) is replaced in the 
expression for tp a by the approximate expression 



kj 



(62) 



2tt 



following from 11 (it when the difference between G^ a 
and Gkj,a is neglected. 
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The expression for the self-energy part 



.4 



Oct 



A 2S 



e_ -B 



B 



5,(63) 



that can be obtained in this case, corresponds to the 
summation of the series 



G = g I ~ goWg I + g Wg Wg I 
where 



(64) 



.9o = 



(uj-e-)- 1 
(u-s A 



W = B 



A a 

A 2S 

-A 2a A 0a 
Ais —Aqo- 



(65) 



A sum of the diagrams with loop-like inclusions into 
the line of the electron single-site Green's function cor- 
responds to this series in a diagrammatic representation. 
Such inclusions lead to the renormalization of energies 
of the electron levels (HE|- In particular, at U = oo 



G^(oj) = ((X^\X"°)) u = J— 

Z7T U 



Aqo 



V4* 

An, 



(66) 



Energy shift 

Ae_ = A 2a B 



A 0a 



t u {Xrxy) (67) 



coincides in this case with the previously obtained one 
in a number of papers (see, for example 15115818 01 using 
a more complicated (in comparison with Hubbard-I ap- 
proximation) decoupling procedure in equations for the 
Green's function G£(w). 

6.3 Hubbard-III approximation (ip a — 0; f2 CT = 

Jo — Ra/Ao a A 2a ) 

We can pass on to this approximation neglecting, at first, 
the renormalization of the atomic electron levels and, 
secondly, approximating 



(X° 



AoaA 2s 



1 



(X 



22 \ 



Aoa-A 2d 



1 



(68) 



in the case of half filling (when (X 00 ) = (X 22 )) in the 
expression that becomes exact only in the U — > 
limit. Consequently, an effective potential of dynamical 
mean field n ff (w) takes the form 



fi CT (w) = J a (w) + Jg(u) — Jg{U — 2fi — oj) 



(69) 



It corresponds (together with the expression l|54Tl for 
the electron self-energy) to the Hubbard-III approxima- 
tion [SS]. A potential includes (besides the coher- 
ent potential J CT (a;)) the terms which describe a scatter- 
ing on the spin and charge fluctuations. Electron energy 
spectrum consists in this case of two subbands only at 
U > U c where the critical value U c corresponds to the 
metal-insulator transition. 

6.4 Alloy-analogy (AA) approximation (R a = 0, 

(p a = 0; 0, a (uj) = JctM) 

The scattering processes are taken here into account only 
via coherent potential. The single-site Green's function 
looks like 



A a 



Ao 



LO — E- — Jo- Ul — S+ — J a 



(70) 



in this case. This expression is analogous to the lo- 
cator function for a binary alloy jHS]- The procedure 
of the G£(w) function calculation corresponds to the 
CPA method. Let us write for this approximation an 
irreducible, according to Dyson, self-energy part T, a = 
uj — E- — SgT 1 , using the expression ll54"jl at Vt a = J a : 



A 2S U 




A aU 



oj + fj, — J a 

Or, after excluding a coherent potential 

A 2S U 



Scr(w) 



l-GS(w)(tf-E a (w)) • 



(71) 



(72) 



This equation corresponds to the alloy-analogy (AA) 
approximation |51| . 

6.5 Modified alloy-analogy (MAA) approxima- 
tion = 0; VL a = J a + V<p a /A 0a A 2 , rT = J a + B) 

An AA-approach is supplemented here by the inclusion 
of renormalization of single-site electron levels. It can 
now be obtained that 




A aU 



(73) 



(74) 



^a=A 2a U , , 

LO — E- — \l 

This relation can be transformed into the equation 

*<»>-W('-%^) 

known in the so-called Modified AA-approach [51.66J. 

One can see from the quoted specific cases that the 
approach developed in this work includes a number of 
known approximations giving in addition their unifica- 
tion and generalization. The proposed scheme is more 
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complete than Hubbard-III approximation (which in its 
turn is the most general of the quoted ones) and dif- 
fers from it by the allowance for a self-consistent renor- 
malization (due to the static internal field) of the local 
energy spectrum as well as by the modification of the 
potential (la constituent parts to a more elaborated in- 
clusion of the magnon and charge fluctuation scattering 
processes. Participation of Bose-particles in such a scat- 
tering is taken into account in our scheme in a more 
consistent way. 

Quantitative changes in the electron spectrum (in par- 
ticular, in the electron density of states) and then in the 
thermodynamics of the model, that might be the conse- 
quence of applying the approach suggested herein, can 
be the subject of subsequent calculations with the use of 
numerical methods. 



7 Simple applications of the method 

Let us demonstrate here the potentialities of the devel- 
oped approximative scheme using the examples of two 
models (the Falicov-Kimball model and the simplified 
pseudospin-electron model) which are analytically solv- 
able in the DMFT approach. 



7.1 Falicov-Kimball model 

The Falicov-Kimball model in its initial version ^2] was 
proposed for the description of the metal-insulator trans- 
formation in compounds with the transition and rare- 
earth elements. The itinerant and localized electrons are 
included into consideration in the model. The simplified 
but sufficiently complete formulation of the FK model 
was given in the set of subsequent publications where 
the Hamiltonian was considered as a specific case of the 
Hamiltonian of the Hubbard model on the assumption 
that the electron transfer from one lattice site to another 
one takes place only in the case of the selected (a =t) 
orientation of spins (see for example |11I90I91I92| 1. Elec- 
trons with the opposite spin orientation (a =J.) remain 
localized. They can effect the energy of delocalized elec- 
trons being a source of scattering. In this case the Hamil- 
tonian of the model is the following [TT1I901I9T] 



E 



tijatt a 3l ■ ( 75 ) 



It should be mentioned that in Q75J1 the itinerant and 
localized particles are of the same nature and possess a 
common chemical potential. Extension of the model on 
the case of the motionless particles of different nature 
(ions, impurities, spins, etc.) results in the Hamiltonian 
of the system of spinless fermions on a lattice which are 
moving in the random field given by the variables (wj = 



ij i 



c,w, 



(76) 



Here the chemical potentials of electrons (/x) and immo- 
bile particles (-E) are different. Among various appli- 
cations the Hamiltonian i|7fit describes the conducting 
electron subsystem of the binary alloy. 

Thermodynamics of the Falicov-Kimball model de- 
scribed by Hamiltonians and Il7fifl was considered 
mainly at the fixed concentration of localized particles 
Pion = jf J2i w i = const (Pion = (^ij.) = const) and in 
the regimes of the given values of the chemical potential 
fi or electron concentration: p e = (rtf) or n = J2 a ( n <?) 

[TTllloll9nll9T] . 

In the DMFT approach the Hamiltonian of the effec- 
tive single-site problem looks in the model ij75ll like 



H eff = (U- 2fi)X 2 



x° 



+ v 

+ H £ . (77) 
In this case, in the equations of motion for A-operators 

[X°\H cS ] = -nX°T + V(A 00 + A T % , 



[X 12 , H cS ] = (U- fj,)X± 2 + V(X 22 + A U )£ T . (78) 

those terms which are responsible for the scattering with 
the participation of Bose-particles (magnons and charge 
excitations) are absent. Only the components 



Z° T =V(A 00 + AiT)£ T 

z 12 = v(x 22 + xii)^ 



(79) 



of irregular parts of the time derivatives of A-operators 
are present. 

The corresponding irreducible Green's functions are 
equal to 



2 7 r((Z T|zT°» L , = 27rT/ 2 A 0T G T M = A 0T J T (w) 
27r((Z i2 \Z 2l )) LU = 2ttV 2 A L2 G 1 (u;) = A i2 J^) 
{{Z^\Z^))^ = {{Z^ 2 \X^)) U = Q. 

Using now the formulae Ij^JI - ftfflfl we obtain 



(80) 



2TiG\ a) {Lo) 



A 



or 



Ul + /! — J-[(u>) 



+ 



A 



12 



ui + \i — IJ — Jf (u>) 



(81) 



This expression has the same structure as the Green's 
function for the AA-approximation and is exact for the 
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Falicov-Kimball model (see, for example, JH]). It can be 
seen directly when we write the equations of motion for 
operators Z ^ , Z^ 2 and take into account two obstacles: 
(i) sums (X 00 + X acr ) and (X 22 + X ss ) are in this case 
an integral of motion and (ii) for the £ CT operator the 
relation 



H eS ] = K a {u>)Z a + V{X°° + X 32 ■ 



(82) 



takes place in the frequency representation; here the 
function K a (w) is connected with the Fourier transform 
of the Green's function ltl9|l by 



(83) 



It should be mentioned that in the auxiliary Fermi- 
field approach one can obtain for the Falikov-Kimball 
model in the framework of the equation of motion 
method an exact expression for the grand canonical po- 
tential fl a of the effective single-site problem. 

Let us represent the exponential operator 
cxp (—(3H c f{) with the Hamiltonian Q77JI in the form 



v = ni e-K Hl+H <) + (1 - Tu)e-* ffa+ff «> , (84) 



where 



Hi — Hi + Hint , H2 — H 2 + Hint , 

fl? = (U- y) n] - n, H\ = —fin-f , 

Hint = V(aU« + iW) . (85) 

By averaging over the £-field and taking the trace over 
the variables <Zf , o| we obtain an operator 

= nM^iP))* + (1 - n 1 )Z O2 (or T 08))§ , (86) 



where erf (/?) is a cr-matrix of the form presented in 1112)) 
with the retarded interaction V 2 Q] (r — r'); (. . .)■' is sta- 
tistical average with the Hamiltonian Hf, 



Z 01 = + e"^-*) 
Z02 = 1 + . 



(87) 



Operator 



-'imp 



-0Hi 



plays the role of the single site ("impurity") statistical 
operator for the electrons with spin |; 



Z imp =^2z 0i (a^(f3))° 



(89) 



For the grand canonical potential we obtain an expres- 
sion 



^imp Q In 2> lm p 



= -Qlnfcr / Z 0i e-™^ , 



(90) 



where the notation j3Qi — — ln(<7f (/?))? is used. 

It is easy to see that the following relations take place 



oft" 



-PQ* 



H2 



"imp 

e /3M + e -/9(l7-/i) 

— 1 

■7. 

-'-'imp 



-PQi 



(91) 



The quantities Q\ and Q2 can be found by means of 
differentiation procedure with respect to the interaction 
constant V. From the one side, using the relations l|90|l 
and Il91fl we get 



dfli 



dV 

From the other side, 



A 



!2 



dQi 
dV 



A 



Of 



dQ 2 
dV 



(92) 



^ E = (^e ff )=2((xto + X^)ef) (93) 

and according to the spectral theorem for Green's func- 
tions 



dQi 



dV 



= 2 



(111.' 



(94) 



The function ((£f |a|)) w is calculated with the help of 
the equation of motion method. Using ll82l) we have 

m\a\)) u =Vg^u)((a 1 \a\)) u; (95) 

and finally, after substitution expression jSH 



Zir 



A 



01 



A 



12 



uj + fi - V 2 Q\ (oj) 



lu + n-u -v 2 gi(uj) 



(96) 



Using this expression in and comparing the ob- 
tained result with I|92J1 we get a result 



+00 



duj 



In 1 - 



J T (w + iO+) 



u> — e+. 



(97) 

It corresponds to the expression which has a form of a 
sum over the Matsubara frequencies 



Q 



1,2 



111,',, 



■ £+.- 



(98) 



and was obtained for the first time in ^T]. Formulae 
(§U for the Green's function and I|90I98|I for the grand 
canonical potential of the single site problem were used 
as a basic expressions for the consideration of the energy 
spectrum and thermodynamics of the FK model in [111 
Ull94H95l in the framework of the DMFT. 
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7.2 Simplified pseudospin— electron model 

In recent years the pseudospin-electron model (PEM) 
has been among the actively investigated models in the 
theory of strongly correlated electron systems. The 
model appeared in connection with the description of the 
anharmonic phenomena in the high-T c superconductors 
and in search of the mechanisms that favor the high val- 
ues of the transition temperatures into superconducting 
state. In addition to the Hubbard type correlation an in- 
teraction with the locally anharmonic lattice vibrations 
(such as vibrations connected with the oxygen sublattice 
ions in the YBa2Cu307_a crystals [97 98 99 ) is included 
into the model. The corresponding degrees of freedom 
are described by the pseudospin variables with S = 1/2. 
The Hamiltonian of the PEM has the form analogous to 
© with 



-.9 



^nurSf-hSf+SlSf (99) 



(see [ED] as wel1 as [TTTnrRMTnSp . Here h is internal 
asymmetry field; SI is a parameter of the tunnelling type 
splitting. 

The model described by the Hamiltonian H99I1 is more 
complicated for consideration than the Hubbard one. 
An analysis of the energy spectrum, thermodynamics 
and charge susceptibility of the model was performed 
in |103I104| in the case U = oo using the generalized ran- 
dom phase approximation (GRPA) |7fill77| . The single- 
electron spectrum in this approach is described in the 
spirit of the Hubbard-I approximation and is splitted 
due to the gnS z interaction into subbands at the any 
value of the coupling constant. 

The DMFT method in its standard formulation based 
on the diagrammatic expansions for the Matsubara 
Green's functions and on the written in the form 111 2ll 
expression for the interaction with the effective field (a 
coherent potential) J a (r — r') was applied by now only 
in the case U = and SI = |H3| - In the GRPA such 
a simplified pseudospin-electron model was considered 
in |105U106j . It is possible to solve analytically an effec- 
tive single site problem for this case too. Let us illustrate 
this with the help of the described above approach. 

The effective Hamitlonian for the simplified model is 
as follows 

H cB = -/xy^ng + n a S z - hS z 

a a 

+ y^ta ff + 4e,) + i { . (ioo) 

a 

Let us write the required single-site Green's function 
in the form 

G^H = «P+a ff |P+at)) w + «p-a ff |p-4)) w , ( 101 ) 

where P ± = 1/2 ± S z are the operators projecting into 
states with a given pseudospin orientation. 



Following the procedure described in section 3 we con- 
sider the equation of motion. 

[^.ffrff] = E ± P ± a a + VPH* , (102) 
where E ± = —fizkg/2. The irregular part in this case is 

z ± = vp±r„ = vp ± s, a . 



The different-time decoupling gives 



((p^\p ± ^)) u = (p ± M <r \C))^ 

As a result, we obtain Dyson equation 

G±(u;) = G± + G±M±Gt 



with 



G 



(P±) 



o 



It follows herefrom that 



(P±) 



Moj) 



(I* 



P+ 



+ 



(P- 



(103) 



(104) 



(105) 



(106) 



(107) 



This expression coincides with the one obtained in the 
d = oo limit in the framework of DMFT [EH- The 
different-time decoupling l|104Jl is also an exact proce- 
dure in this case. 

Let us mention that the two-pole structure of the 
single-site Green's function leads to the effect of the 
metal-insulator transition type at t ~ g. In the case 
t > g the electron spectrum consists of one broad band 
while at t < g there appears a gap and the splitting 
into two Hubbard-type bands takes place |BU ■ It should 
be stressed on this occasion that the Hamiltonian of the 
simplified PEM corresponds to the Hamiltonian jS§ of 
the FK model at the replacement Sf — Wi — \. There 
exists however an essential difference in the regimes of 
thermodynamical averaging: in the PEM the value of 
the field h is fixed (which is an analogue of the chem- 
ical potential — E in l|76Jl L but not the mean value of 
pseudospin. The special features of the electron and 
pseudospin subsystems behavior and phase transitions 
in PEM in this case are investigated and described in 

|S4ll1oMI1ofi| . 

8 Perturbation theory in terms of elec- 
tron hopping 

Presented in the previous sections an irreducible Green's 
function approach allows to construct various approxi- 
mations for the single-electron properties, but does not 
allow in general case to describe the thermodynamics 
in a self-consistent way. In this section we present a 
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different approach based on the Wick's theorem and di- 
agrammatic technique for the Hubbard operators |7()j . 

We consider the lattice electronic system that can be 
described by the following generalized statistical opera- 
tor: 







a(/3) = Texp I- Jdr J At' 
I o o 



where 



Hq — 22 Hi 



(109) 



is a sum of the single-site contributions and for the Hub- 
bard model we have 

Hi = Un^riii - p{n^ + n a ) - h(n^ - n^) , 

t^T-T^^t^T-T') . (110) 

In addition, for the Falicov-Kimball model 



_ , U 3 S{t-t') forcr=T 

** K ' i fora=| 



(111) 



On the other hand, for the auxiliary single-site problem 
(JEl of the DMFT we must put 



=5 lJ J a (T-T r ) 



(112) 



It is supposed that we know eigenvalues and eigen- 
states of the zero-order Hamiltonian l|109|) . 



Hi\i,p) = X P \i,p) 



(113) 



and one can introduce Hubbard operators in terms of 
which zero-order Hamiltonian is diagonal 



(114) 



For the Hubbard model we have four states \i,p) = 
\i,ritf,riii) J^njl: |i, 0) = |i, 0, 0) (empty site), \i, 2) 
/. 1.1 (double occupied site), \i, t) = \i, 1,0) and 
) = \i, 0, 1) (sites with spin-up and spin-down electrons) 
with energies 



A = , X 2 = U-2p, 
Aj. = h - p , Af = —h — p 



(115) 



The connection between the electron operators and the 
Hubbard operators is the following: 



Xf + Xf CT ; at* = X° a + aXf 



(116) 



Our aim is to calculate the grand canonical potential 
functional 



Q = -imSpp = n -^m(<7(/?)>0, 

O = -ilnSpe-^° , 

P 

single-electron Green's functions 
G l3a (T t') = (Tal{T)a ja {T')) = ^, ( f-T ') 



(117) 



7T (H8) 



and mean values 



1 dfl 



N dp a 1 



(119) 



where p a = p. + oh is a chemical potential for the elec- 
trons with spin a. Here, (...) = Sp(. . .p)/Z, Z = Spp, 
or in the interaction representation 



1 



-<...a(/3))o = (...a(/3))oc, (120) 



V "' (HP))' 

where (. . .) = Sp(. . . e-^ H °)/Z ; Z = Spe"^ . 

We expand the scattering matrix <r(/3) in 111081) into 
the series in terms of electron hopping and for (er(/3))o we 
obtain a series of terms that are products of the hopping 
integrals and averages of the electron creation and an- 
nihilation operators or, using lllfijl . Hubbard operators 
that will be calculated with the use of the corresponding 
Wick's theorem. 

Wick's theorem for Hubbard operators was formulated 
in Oj (see also Ref. [J^l and references therein). For 
the Hubbard model we can define four diagonal Hub- 
bard operators X pp (p = 0, 2, J., f ) which are of bosonic 
type, four annihilation X oi , X°\ X^ 2 , X 12 and four 
conjugated creation fermionic operators, and two anni- 
hilation X^ , X 02 and two conjugated creation bosonic 
operators. The algebra of X operators is defined by the 
multiplication rule 



the conserving condition 



and the commutation relations 



[XT>,X?] ± = 6 ii (8 v X?±6r q X?) 



(121) 



(122) 



(123) 



where one must use anticommutator when both oper- 
ators are of the fermionic type and commutator in all 
other cases. So, commutator or anticommutator of two 
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Hubbard operators is not a c number but a new Hub- 
bard operator. Then the average of a T products of X 
operators can be evaluated by the consecutive pairing, 
while taking into account standard permutation rules 
for bosonic and fermionic operators, of all off-diagonal 
Hubbard operators X pq according to the rule (Wick's 
theorem) 



xr*( Tl )X™ (r)=-*oiWr-n) 

xWinlXFWh (124) 

until we get the product of the diagonal Hubbard op- 
erators only. Here we introduce the zero-order Green's 
function 



9pq 



V 

- p (t-ti)a p , (±n±(X pq ) 



T > Tl 



±n±(\ pq ) - 1 T < Tl 



(125) 



where X pq = \ p — X q and n±(X) = (c^ A ± l) , and its 
Fourier transform is equal 



ko v — X pq 



(126) 



Applying such pairing procedure to the expansion of 
(a(/3))o we get the following diagrammatic representa- 
tion: 



exp-! 



-(.127) 



-00-cfe--CX>- 

where arrows denote the zero-order Green's functions 
ijl2fi|l . wavy lines denote hopping integrals and □, ... 
stay for some complicated "n vertices", which for such 
type perturbation expansion are an irreducible many- 
particle single-site Green's functions calculated with the 
single-site Hamiltonian 11141 . Each vertex (Green's 
function) is multiplied by a diagonal Hubbard opera- 
tor denoted by a circle and one gets an expression with 
averages of the products of diagonal Hubbard operators. 

For the Falicov-Kimball model expression l|127H re- 
duces and contains only single loop contributions 



where 



exp< 




-©- 



P^~ = hii; P = 1 — hii, fi* = fi — U/2 and by intro- 
ducing pseudospin variables Sf = (P 4 — P i _ )/2 one can 



transform the Falicov-Kimball model into an Ising-type 
model with the effective multisite retarded pseudospin 
interactions. Expression 112811 can be obtained from the 
statistical operator Q108J1 by performing partial averag- 
ing over fermionic variables, which gives an effective sta- 
tistical operator for pseudospins (ions). 

So, after applying Wick's theorem our problem splits 
into two problems: (i) calculation of the irreducible 
many-particle Green's functions (vertices) in order to 
construct expression 112711 and (ii) calculation of the av- 
erages of the products of diagonal Hubbard operators 
and summing up the resulting series. 



9 Irreducible 
functions 



many-particle Green's 



For the Hubbard model by applying the Wick's theorem 
for X operators one gets for two-vertex 



= g a0 (u u )(Xr + X?°) 

+ 92a {U V ) (X? 2 + Xf & ) 



(129) 



for four-vertex 



U) v (7 



A (4) - v 

1(7(7 V 



^>v+mi ^v'+m, ^V) (130) 

= A > 4 00 g <TO (w 1 ,)g CT o(w 1 ,+m) (U + U 2 g 2 o(uv+v>+ m )) 

X<foo(<*v)<feo(<*V+m) 

+X 22 g 2s {uJu)g2s{^u+m) (U - U 2 g 2 o{u v + v >+ m )) 

X.g2a{Uv')g2o{Uv>+m) 

+X[° ' g a o{^u)gcxo{^u+m) {U + U 2 g< T& (u v - v >)) 

Xg2a{u u ')g2a{Uv'+m) 

+X°° j2sK)S25K+m) (U — U 2 g a a{u v -v')) 

xg&o{u U ')gvo{uv'+ m ) , 



A 



(4) 







and so on. Expressions l|129|) and H130|l and for the ver- 
tices of higher order possess one significant feature |70 ; . 
They decompose into four terms with different diago- 
nal Hubbard operators X pp , which project our single- 
site problem on certain "vacuum" states (subspaces), and 
zero-order Green's functions, which describe all possible 
excitations and scattering processes around these "vac- 
uum" states: i.e., creation and annihilation of single elec- 
trons and of the doublon (pair of electrons with opposite 
spins) for subspaces p — and p — 2 and creation and 
annihilation of single electrons with appropriate spin ori- 
entation and of the magnon (spin flip) for subspaces p =f 
and p =1. 

In compact form expressions l|129ll and 113011 can be 
written as 



9(7(p){Uv) 



(131) 
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and 

]0[ = ^ Xi P 9a(p) (uv)g<j( P ) {u y+m )U aS(p) (uv, u v > \ui m ) 



X 9a(p){^v')gs(p){^v' +m) 



where 



g a o{LOv) for p = 0, a 
92a{u v ) for p = a,2 



(132) 



(133) 



Here 



^ff(p)( w l/,W„'|w m ) 

= (U ± t/ 2 .g 2 o(^+^+m) for p = 0, 2 
\J7 ± U 2 g aS (tjJv-v* ) for p = a, a 

Uaa(p){u Vl LOv>\LO m ) — Uv<r(p) (oV , UJ V \ L0 m ) 



(134) 



is a renormalized Coulomb interaction in the subspaces. 
In diagrammatic notations expressions H130J1 or I|132J1 

can be represented as 



n X 



± 



3 2 



for p = 0, 2 



for p ■■ 



(135) 

where dots denote Coulomb correlation energy U = A 2 + 
A — A | — A| and dashed arrows denote bosonic zero-order 
Green's functions: doublon <?2o( w m) or magnon g a a{^m)- 
Expression for six-vertex contains the contributions 
which can be presented by the following diagrams: 



(136) 



with the internal vertices of the same type as in Q135J1 
and contribution which can be presented diagrammati- 
cally as 

(137) 



>V-< 



So, we can introduce primitive vertices 



x>- X 



(138) 



by which one can construct all n vertices in expansion 
(I127H according to the following rules: 

1. n vertices are constructed by the diagonal Hubbard 
operator X pp and zero-order fermionic and bosonic 
lines connected by primitive vertices H138H specific 
for each subspace p. 

2. External lines of n vertices must be of the fermionic 
type. 



3. Diagrams with the loops formed by zero-order 
fermionic and bosonic Green's functions are not al- 
lowed because they are already included into the 

formalism, e.g., gives r~~^. 

For n vertices of higher order a new primitive vertices 
can appear but we do not check this due to the rapid in- 
crease of the algebraic calculations with the increase of n. 
Diagrams H135J1 , HI 3611 , and H137H topologically are trun- 
cated Bethe-lattices constructed by the primitive ver- 
tices i|138H and can be treated as some generalization of 
the Hubbard stars |l()7ll()8ll()9j in the thermodynamical 
perturbation theory. 

It should be noted that each n vertex contains 
Coulomb interaction U as in primitive vertices I|138J1 (de- 
noted by dots) as in the denominators of the zero-order 
Green's functions II126II . In the U — * oo limit, each term 
in the expressions for n vertices can diverge but total 
vertex possesses finite U — > oo limit when diagrammatic 
series of Ref. [23 are reproduced. 

The second problem of calculation of the averages of 
diagonal X operators is more complicated. One of the 
ways to solve it is to use semi-invariant (cumulant) ex- 
pansions as was done in Refs. [211 an d EH f° r the U = oo 
limit. Another way is to consider the d — oo limit where 
new simplifications appear. 



10 Dynamical mean-field theory 

In general, the grand canonical potential for lattice is 
connected with the one for the auxiliary single-site prob- 
lem of the DMFT by the expression pi] 



lX>G,K,k)j . 



(139) 



On the other hand, we can write for the grand canonical 
potential for atomic limit fli mp the same expansion as 
in H127J1 but now we have averages of the products of 
diagonal X operators at the same site. According to 
l|121l) we can multiply them and reduce their product 
to a single X operator that can be taken outside of the 
brackets and exponent in I|127J1 and its average is equal 
to 



(X pp ) 



E e -/JA, 



Finally, for the grand canonical potential in atomic limit 
we get 



5 -/3f2 (p) 



(140) 
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where 



0( p ) = X p + 



P 



are the "grand canonical potentials" for the subspaces. 

Now we can find single-electron Green's function for 
atomic limit 

G W( r _ r ') =/3 m ™*> 



8J a {T-T') 

=Y, w p g °(p)( t - t 1 ) ' ( 142 ) 



where 



G CT ( p )(r-r') = /?- 



sn 



■*Wr-^) <I43) 
are single-electron Green's functions for the subspaces 
characterized by the "statistical weights" 

e -/30 (p) 

™* - ( 144 ) 

and our single-site atomic problem exactly splits into 
four subspaces p = 0, 2, J., f . 

We can introduce irreducible parts of Green's func- 
tions in subspaces S CT (p) (^v) by 



1 



where 



(145) 



(146) 



According to the rules of the introduced diagrammatic 
technique, n vertices are terminated by the fermionic 
Green's functions [see H135I) . \ 13611 . and (113 T|> ] and this 
allows us to write a Dyson equation for the irreducible 
parts and to introduce a self-energy in subspaces 



(147) 



where self-energy Yi a ^{oj v ) depends on the hopping in- 
tegral J a i (uv ) only through quantities 

1 3V(p)(<<V) = Ga'{p)(^') - S CT '(p)(<<V) 

= 3 ff'( P )(^') J <T'(^') 

x {1 + E^KOJ^av) + • ■ ■} ■ (148) 

It should be noted, that the total self-energy of the 
atomic problem is connected with the total irreducible 
part by the expression 



(149) 



and it has no direct connection with the self-energies in 
the subspaces. 

The fermionic zero-order Green's function 113311 can 
be also represented in the following form 



where 



9a(p) 



l (0) 



1 



for p = 0, a 

1 for p = 2, a 



(150) 



(151) 



is an occupation of the state \p) by the electron with spin 
a, and Green's function I|145J1 can be written as 



G<r(j>) (u v ) = \u v + fi a - Uny 



(/--) 



(152) 



Now, one can reconstruct expression for the grand 
canonical potentials fi/ p ) in subspaces from the known 
structure of Green's functions. To do this, we scale hop- 
ping integral 

J a {u v ) -» aJ„(w„) , «e[0,l], (153) 

which allows to define the grand canonical potential as 



(p) — 



J da jj^2 Ja(^u)G a{p) (uj v ,a) (154) 



and after some transformations one can get 

fi(p) = A p - - In [1 - J (T (^)S (T(p) (w„j\ 

3E^)W*'(p)W +i (p)' ( 155 ) 



p 



where 



*(p) = 



1 



(156) 



is some functional, such that its functional derivative 
with respect to \? produces self-energy: 



P 



5$ 



(p) 



(157) 



So, if one can find or construct self-energy S .( p )(w l/ ) he 
can find Green's functions and grand canonical poten- 
tials for subspaces and, according to 1114011 and 1142L 
solve auxiliary single-site problem. 
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Starting from the grand canonical potential l|14fljl and 
11551 one can get for mean values 1119|) . 



,(°) 



\ [ G ^(p)(^) - -^(p)M] 



(p) 



(158) 



where in the last term the partial derivative is taken over 
the \x a not in the fields * cr (p)(aj 1/ ) l|148|) . The second 
term in the right-hand side of l|158|) can be represented 
diagrammatically as 



and the first contributions into the last term are follow- 
ing 




(160) 

where double lines denote quantities ^> a ^ p j(uj v ). Loop 
is connected with the superconducting or magnon 
susceptibilities for subspaces p — 0, 2 or p — a, a, respec- 
tively. 

For the single atom = 0] we have = 0, 

Go{p){u v ) = 3 (T ( p )(a; 1/ ) = g a ^{u} u ), and 



(0) 

w pK(p) > 



(161) 



but in the general case [J CT (a-v) ^ 0] we cannot prove 
that the sum rule 



(162) 



is fulfilled. 



71 f 



n l =w 2 +w l , 



(166) 



which immediately gives results of ( see a ^ so 

Ref. jHU)- Expression I|165J1 is the Matsubara represen- 
tation of the Green's function 181J1 and the second term 
in l|164JI corresponds to the quantities Qi, 2 l|97ll . which 
give its real axis representation. 

For the Hubbard model there are no exact expression 
for self-energy but the set of Eqs. 114511 . 114711 . and 1155jl 
allows one to construct different self-consistent approx- 
imations. 

10.2 Alloy-analogy approximation 

The simplest approximation, which can be done, is to 
put 

S CT(p) K) = 0; $(p)=0 (167) 



which gives 



(168) 



and 



fi(p) = A p - — In [l - J a (uj v )g^ p )(uj v )] (169) 



and for the Green's function for the atomic problem one 
can obtain a two-pole expression 



Wq + W a 



1U) V — Ao-0 — Ju\Wv) 

W 2 + W 9 
kd„ — \ 2S — J<t(ll>is) 



(170) 



of the alloy-analogy solution for the Hubbard model, 
which is a zero-order approximation within the consid- 
ered approach. For this approximation, mean values 
H119I) are equal to 



10.1 Falicov— Kimball model 

For the Falicov-Kimball model Ji(oj v ) — and accord- 
ing to tup 



s t(p)( w ") = > s t(p)( w ") - 5t(p)( w ^) ; 



$ (p) ee0 



and 



G^K): 



1-71 



\u v - A T0 - J]{u„) 

™L 

\u v - X 2l - J](u v ) 



(163) 



fi(p) =A p --^hi[l- J-[(u v )g^(u v )] , (164) 



(165) 



1 r<(°>) ( \ A , w + w a w 2 +w s 

n a = -g > J G^> {w v ) + w 2 + w a - 



(171) 



and, for some values of the chemical potential, they can 
get unphysical values: negative or greater then one. 

10.3 Hartree Fock approximation 

The next possible approximation is to take into account 
the contribution from diagram I|159J1 and to construct 
the equation for the self-energy in the following form: 

Ea(p)K) = ^5^tr*ff( p )(w^) ) (172) 
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which, together with the expression for mean values 



(o) 



r(v) - ,l a(p) f p 



,(0) 



1 1 







(p) - 2 ^ 2 tanh 2 ^ nff (p) ~ 



(173) 



gives for the Green's function in the subspaces expression 
in the Hartree-Fock approximation: 



1 



+ fJ.tr- Un s{p ) - J a {uj u ) 



(174) 



Now, grand canonical potentials in the subspaces are 
equal 

fi(p)=A p - 4^I ln t 1 - ^(^) s <t(p)(^)] 

1/(7 

-U (r Mp) - (n ff(j>) - 4°^)) . ( 175 ) 



$ (p) = U[n 



l c{p) 



1 Up)) ( 



(0) 

n <y{p) n s(p) 



and for the single-site Green's function 114211 one can 
obtain a four-pole structure 



V + fl<r - Un<r(p) - J a (0Jv) 



(176) 



Expression II176I) . in contrast to the alloy-analogy solu- 
tion I|17()J1 . possesses the correct Hartree-Fock limit for 
small Coulomb interaction [/<(: 



1 



\uj v + [L a — Un s — J a {u v ) 



(177) 



when w p « 1/4 and n a ^ w n CT = 6 G^ On 
the other hand, in the same way as an alloy-analogy 
solution, it describes the metal-insulator transition with 
the change of U. 

In Fig. the frequency distribution of the total spec- 
tral weight function 



(178) 



as well as contributions into it from the subspaces [sep- 
arate terms in 1117611 ] are presented for the different 
electron concentration (chemical potential) values. One 
can see, that the spectral weight function contains two 
peaks, which correspond to the two Hubbard bands. 
Each band is formed by the two close peaks: p = 
and a for the lower Hubbard band and p = 2 and a for 
the upper one, with weights w p <I144|) . The main con- 
tributions come (see Fig. |2J from the subspaces p = 
for the low electron concentrations (n < 2/3, (j < 0), 



p = 2 for the low hole concentrations (2 — n < 2/3, 
ji > U) and p — a, a for the intermediate values. For 
the small electron or hole concentrations, the Green's 
function for the atomic problem 117611 possesses correct 
Hartree-Fock limits too. 

Such four-pole structure of the single-electron Green's 
function can be obtained also for the one-dimensional 
chain with the N — 2 periodic boundary condition (see 
Appendix in Ref. |70|). which is equivalent to the two- 
site problem considered by Harris and Lange |58j . Here, 
two poles correspond to the noninteracting electrons or 
holes, which hope over the empty sites, and give the 
main contribution for small concentrations. The other 
two poles give the main contribution close to half-filling 
and correspond to the hopping of the strongly-correlated 
electrons over the resonating valence bond (RVB) states. 

So, one can suppose that the Hubbard model describes 
strongly-correlated electronic systems that contain four 
components (subspaces). Subspaces p = and p = 2 
describe the Fermi-liquid component (electron and hole, 
respectively) which is dominant for the small electron 
and hole concentrations, when the chemical potential is 
close to the bottom of the lower band and top of the 
upper one. On the other hand, subspaces p =t and J. 
describe the non-Fermi- liquid (strongly correlated, e.g., 
RVB) component, which is dominant close to half- filling. 
The plateau at half filling for w p (p =T,I) can be as- 
sociated with the antiferromagnetic phase. Within the 
considered Hartree-Fock approximation, at n ps 2/3 
and 2 — n ~ 2/3, we have transition between these 
two regimes: Fermi liquid and non-Fermi liquid. It 
reminds us the known properties of the high-T c com- 
pounds, where for the nondoped case (n = 1) com- 
pounds are in the antiferromagnetic dielectric state, then 
for small doping the non-Fermi-liquid behavior is ob- 
served (underdoped case n < 1) and after some optimal 
doping value, the properties of the compound sharply 
change from the non-Fermi to the Fermi liquid (over- 
doped case). 

The results presented in Figs. and are obtained 
for relatively high temperature. With the temperature 
decrease, on the one hand, the transition between the 
Fermi and non-Fermi liquid becomes sharp and, on the 
other hand, for some chemical potential values there can 
be three solutions of H173JI with two of them correspond- 
ing to the phase-separated states. The consideration of 
the phase separation in the Hubbard model is not a topic 
of this paper and will be the subject of further investi- 
gations. 

At low temperatures, besides the plateau on the con- 
centration dependence of w p for p = a, a at half filling, 
also the plateau for the statistical weights of subspaces 
p = 0, 2 are developed at low electron and hole con- 
centrations, see Fig. 03 The p = and p — 2 compo- 
nents for the low electron or hole concentrations are in 
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the ferromagnetic state, while the non-Fermi-liquid one 
is antiferromagnetic (AF) close to half- filling jlldj . For 
the intermediate concentration values the picture is very 
complicated, even frustrated. It is due to the fact that 
equations for the mean values Ijl73|) have several solu- 
tions in this region, which, on the other hand, are mu- 
tually connected with the dynamical mean field J a (uj v ). 
It is difficult to determine the ground state for this, pos- 
sibly "pseudo-gap", region, which is located between the 
ferromagnetic and antiferromagnetic phases. 

In Fig. 31 we presented the phase diagram (T, U) - the 
temperature of the AF ordering vs correlation energy 
U, which is in a qualitative agreement with the results 
of Refs. |21Ulllll27i] and reproduces the results of the 
Hartree-Fock theory and mean field approximation for 
U <C t and U ^> t, respectively. Our results for the 
AF critical temperature for small U are higher then the 
one of the Quantum Monte Carlo simulations j2l] by 
about a factor at three that describes the reduction of 
the Hartree-Fock solution by the lowest order quantum 
fluctuations |lf)7| . 

10.4 Beyond the Hartree-Fock approximation 

Self-energy in the Hartree-Fock approximation [see 
Eq. 117411 ] describes some self-consistent shift of the ini- 
tial energy levels and does not depend on the frequency. 
All other improvements of the expression for self-energy 
add the frequency dependent contributions. To see this, 
let us consider the contribution into the mean values 
from the first diagram in IjlfiOJI . This diagram originates 
from the following skeletal diagram 



(179) 



in the diagrammatic expansion for functional &( p )- On 
the other hand, such a skeletal diagram produces addi- 
tional contribution into the self-energy 

which is frequency dependent. Also, in order to get a 
self-consistent set of equations, we introduce renormal- 
ized bosonic Green's functions 




D20(p) (^m) 



Arff(p) i^rn) 



1 



\uo m - A 



20(p) 



i^m — ^aa(p) 



^20(p) — ^20 + U— *cr(p)(^) ) 
va 



(181) 



Finally, for the Green's function I|152J1 we get the gen- 
eral representation 



G<j(p){uv) = iu v + [i a - Un s{p) 



(182) 



where the Hartree-Fock contribution f7n S ( p ) is extracted 
and S(T(p) is a frequency dependent part of the self- 
energy, which within the considered approximation is 
equal 

£<7(p)(k>i/) = 

U 2 

±^^D aS(p) {uj u ,Lo ul )^ g^ujv,) , (183) 



where 
D 



Dia{p){uu+v>) for p= 0,2 
D aS ( p) (u v - v > ) for p = a, a 



(184) 



Now, mean values 115811 are equal 



,(0) 



Hp) t p 



(185) 



and for the functional l|15fijl in the grand canonical po- 
tentials in the subspaces we obtain the following expres- 
sion 



vv' 

x *a(p)(^)*<T(p)(kv) • 



(186) 



Expression l|18fi|l besides the Hartree-Fock contribution 
HI 7511 contains also the contribution from the skeletal 
diagram (11791) . 

In order to analyze the structure of the poles in 
l|182ll . an analytical continuation of the expression for 
E CT ( p )(a;„) from the imaginary axis to the real one should 
be done. To do it, we use the well-known identity 



-Y- 



icj, y — A 



= ±n±(A) , 



(187) 



which follows from 112511 . and analytical properties of 
the Green's function 



GJz) 



dw 



3G CT (w-i0+) 



z — U) 



(188) 
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Green's functions in the subspaces G CT ( p )(z), irreducible 
parts S CT ( p )(z), and dynamical mean-field J a (z) all pos- 
sess the same analytical properties. Finally, we get the 
following expressions: 



+ OC 



E a(p )(z) = ± — dujn + (u) ^ = 

71" J Z + UJ- A 2 0(p) 

— OO 

T«-(A 20 ( P )) ^ 2 *<t(p)(A 2 o( p ) - z) (189) 
for subspaces p = 0, 2 and 

W) = ±- Aw.) ^)(- i0+ ) (190) 

n J_ Z-W- Xaa(p) 

r2, 



=F 



for p = o", tr. Analytical continuation of expressions 
118511 and i|186ll can be done in the same way. One can 
see, that contributions H189J1 and l|19fljl diverge in the 
paramagnetic phase close to half filling when A 20 ( p ) = 
and \ a s(p) — 0, respectively, which is an unphysical re- 
sult. 

So, we cannot include into the consideration only one 
contribution from diagram i|179Jl but one have to con- 
sider, besides the fermionic loops, also the bosonic ones 
|112j which correspond to the creation and annihilation 
of the doublons (pairs of electrons), described by the X 20 
and X 02 operators, for subspaces p — 0, 2 and magnons, 
described by the and X^ operators, for p =t,|- 
The such loop contributions of bosonic excitations can 
be summed up and one can obtain 



x (l ! I'D, 



U 



where 



Xaa(p) (Urn) 



, (191) 

(192) 



£ M^&ip) 



for subspaces p = 0, 2 and 

D aS { p ) (Wm) = Daa(p) ( w m) , 



(193) 



(?) 



for subspaces p = <j,<t. Expression 111861) is the first 
term of the expansion of functional H191J1 in the series 
over Xos{p){u m )- 

Now we obtain for mean values the following expres- 
sion 



Mp) 



<r(p) 



£*-(?) 



(194) 



^ 1 U 2 Dl- {p) (uj m )x aS{p) (u; m ) 

(Wm) 



and self-energy contains the frequency dependent part 

= C 7 («ff(p) ~ n f(p)) + £<r(p)M > (19S) 



1 ± J7L' CTff (p)(w m )J Xo-ff( P )(w m ) ± Dto(p)K) 



l-tf(l±tf£W( P )(w (w m ) 



*a(p)K-i-)i for P = 0,2 

*<t(p) (w„_ m ), for p = cr, a ' 



(196) 



that describes the contributions from the doublons 
(charge fluctuations) for the Fermi liquid component 
(p = 0, 2) and magnons (spin fluctuations) for the non- 
Fermi liquid one (p =|, |) with the renormalized spec- 
trum determined by the zeros of denominator in j 19611 . 

Expression Ill91fl for functional has the same form 
as the correction to free energy in the theory of the self- 
consistent renormalization (SCR) of spin fluctuations by 
Moriya |113j . But in our case it describes contributions 
from the single-site bosonic (spin or charge) fluctuations 
with specific renormalization functions different for dif- 
ferent subspaces. Spin fluctuations give the main contri- 
bution close to half filling in the non-Fermi liquid regime 
but for small electron (n <C 1) or hole (2— n <C 1) concen- 
trations the contributions from the charge fluctuations 
must be taken into account. 

11 Concluding remarks 

An analytical approaches for the solution of the effec- 
tive single site problem in the DMFT method for the 
Hubbard-type models described in this article are based 
on the strong coupling scheme that considers the strong 
local interaction as reference system. For the first one, 
it corresponds to the selection of the Hubbard operators 
as basis for the projection procedure for Green's func- 
tions while in the second one the perturbation theory 
over electron hopping is used. Both of them have their 
advantages. 

The equation of motion method together with the av- 
eraging over the auxiliary Fermi field gives an approxi- 
mate interpolating scheme that in specific cases includes 
a number of known approximations for the Hubbard and 
similar models. An examples where the proposed ap- 
proach gives exact results are given (Falicov-Kimball 
and simplified pseudospin-electron models). 

At the same time, the applied procedure of the irre- 
ducible Green's functions introduction and different time 
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decoupling appears to be too simple to obtain the 4-pole 
structure for the single-site electron Green's function. 
An inclusion only of the Fermi-type single-site Hubbard 
operators in the basis at the formulation of the equations 
of motion produces the 2-pole Green's function and only 
the extension of the basis and application of the pro- 
jection and decoupling procedures to the higher order 
functions probably can be able to reveal the more com- 
plicated structure of function Ga\u). Besides, this way 
requires the consideration of the retarded effective inter- 
actions formed by the auxiliary £-field. 

Nevertheless, it should be mentioned, that the sim- 
plicity and accessibility of such approach based on the 
equation of motion method makes it attractive for the 
approximate analytical considerations. It seems useful 
to apply it to the problems which have been considered 
up to now by means of numerical methods (or can be 
solved exactly only numerically). It should be noted, 
that the calculation of the electron mean occupation val- 
ues (and derivation of the equation for the chemical po- 
tential) , as well as the determination of the grand canon- 
ical potential within the equation of motion scheme for 
the Green's functions are elucidated only partially in this 
work. It will be the subject of a separate publication. 

The second approach considered in this article uses 
for the Hubbard-type models a finite-temperature per- 
turbation theory scheme in terms of electron hopping, 
which is based on the Wick's theorem for Hubbard op- 
erators and is valid for arbitrary values of U (U < oo). 
Diagrammatic series contain single-site vertices, which 
are irreducible many-particle Green's functions for un- 
perturbated single-site Hamiltonian, connected by hop- 
ping lines. Applying the Wick's theorem for Hubbard 
operators has allowed us to calculate these vertices and 
it is shown that for each vertex the problem splits into 
subspaces with "vacuum states" determined by the diag- 
onal (projection) operators and only excitations around 
these "vacuum states" are allowed. The vertices possess 
a finite U — > oo limit when diagrammatic series of the 
strong-coupling approach are reproduced. The 

rules to construct diagrams by the primitive vertices are 
proposed. 

In the limit of infinite spatial dimensions the total aux- 
iliary single-site problem exactly (naturally) splits into 
subspaces (four for Hubbard model) and a considered an- 
alytical scheme allows to build a self-consistent Baym- 
Kadanoff-type theory for the Hubbard model. Some an- 
alytical results are given for simple approximations: an 
alloy-analogy approximation, when two-pole structure 
for Green's function is obtained, which is exact for the 
Falicov-Kimball model, and the Hartree-Fock-type ap- 
proximation, which results in the four-pole structure for 
the Green's function. Expanding beyond the Hartree- 
Fock approximation calls for the considering of the fre- 
quency dependent contributions into the self-energy con- 



nected with the self-consistently renormalized spin and 
charge fluctuations. 

In general, the expression 



(197) 

gives an exact four-pole structure for the single-site 
Green's function of the effective atomic problem. In 
H147H zero-order Green's functions 11330 are the same 
for the subspaces p = 0, a and p = 2, a, respectively, 
and correspond to the two-pole solution of the one-site 
problem without hopping. Switching on of the electron 
hopping splits these two poles and the value of splitting 
is determined by the values of the self-energy parts in 
the subspaces, which describe the contributions from the 
different scattering processes. Alloy-analogy approxima- 
tion neglects such scattering processes (E ff ( p )(w„) = 0) 
which results in the two-pole structure for the Green's 
functions 11701 . But, in general, Green's functions pos- 
sess four-pole structure and even the Hartree-Fock ap- 
proximation H176J1 clearly shows it. 

It should be noted that the four-pole structure of the 
Green's function for the atomic problem might not re- 
sult in the four bands of the spectral weight function (see 
Fig. nj. The presented consideration allows us to sup- 
pose that each pole describes contributions from the dif- 
ferent components (subspaces) of the electronic system: 
Fermi liquid (subspaces p = 0, 2) and non-Fermi liquid 
ip =T> I)) an d for small electron and hole concentrations 
(n < 2/3 and 2 — n < 2/3) the Fermi-liquid compo- 
nent gives the main contribution ("overdoped regime" of 
high-T c 's), whereas in other cases the non-Fermi liquid 
one ("underdoped regime"). 
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Figure 1: Spectral weight function p a (uj) 1I178J1 : total and for each subspace, for the different chemical potential 
values: (a) fi = U/2, n = 1; (b) // = -1, n = 0.07; (c) fi = 0.01, n = 0.72; (d) fi = -0.01, n = 0.66 (U = 4, T = 0.2) 
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Figure 2: Statistical weights of the subspaces w p II 14411 as functions of the electron concentration (U = 4, T = 0.2) 
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Figure 3: Statistical weights of subspaces w p and ferromagnetic rap and antiferromagnetic tuaf order parameters 
vs electron concentration for U = 1.56, T = 0.14 
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Figure 4: Phase diagram (T, U) at half-filling n = 1 (AF - antiferromagnetic phase, PM - paramagnetic phase) 
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